Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences

Early studies of synaptic transmission by Bernard Katz and colleagues suggested that neurotransmitter release at graded-potential synapses occurs through statistically independent (i.e. Poissonian) quanta [1, 2]. Subsequent experimental work supported this framework [3]. However, these measurements were performed in vitro on relatively simple synapses and under non-physiological conditions, often converting spiking neurons into graded-potential neurons through the use of channel blockers. Relying on the conventional assumption that vesicle exocytosis follows a Poisson process, measurements of the contrast power transfer spectrum and noise power spectral density of large monopolar cells (LMCs) in the blowfly C. vicina imply a sustained vesicle release rate exceeding 105 vesicles per second per LMC. Given the physical dimensions of photoreceptors and synaptic vesicles, such a release rate appears physiologically implausible. If vesicle release is more temporally structured, low-frequency noise could be suppressed, substantially reducing the vesicle release rate required to account for experimental observations. The reduction of noise at low frequencies is especially advantageous given inputs such as photoreceptor signals which are already low-pass filtered. Visual activity generates substantial extracellular potentials within the lamina cartridge [4]. We propose that these extracellular potentials regulate vesicle release by modulating the voltage sensors that trigger exocytosis. We provide experimental evidence for the connection between currents driving the LMC and the extracellular potentials during visual activity, and demonstrate, using simple models, how effective "Poisson" rates are maximized due to vesicle regularization.

Surface electromyography (sEMG) is the most practical non-invasive interface for myoelectric prostheses, exoskeletons, and rehabilitation systems, but power-line interference (PLI) contamination and excessive digital pipeline group delay still limit its clinical adoption. This paper proposes a co-designed analog-digital correction system combining a high-CMRR front-end with an exponentially-windowed RMS (EMRMS) envelope estimator and a recursive single-tone PLI canceller. We present a closed-form CMRR model capturing the electrode-skin imbalance, and provide a complete stability analysis of the LMS canceller. The EMRMS estimator reduces the computational overhead from[O] (L) to strictly[O] (1) in both time and space complexities. Featuring no data-dependent branching, the algorithm achieves deterministic algorithmic execution time (zero jitter under an RTOS environment) and is natively compatible with fixed-point arithmetic on microcontrollers lacking a hardware Floating-Point Unit (FPU). A reference implementation reaches an 8.2 {micro}s median per-sample latency, yielding an end-to-end delay of[~] 30 ms -- leaving a generous >90 ms budget for electromechanical actuation -- while requiring an active CPU duty cycle of merely 1.6%, enabling prolonged deep-sleep intervals. Validation on the public Ninapro DB2 dataset demonstrates a 13.9 dB mean SNR improvement (averaged across 12 channels; single-channel comparison: 9.7 dB, Table 3) and a 70.0 {micro}V envelope RMSE against a length-200 rectangular reference. Paired Wilcoxon signed-rank tests confirm statistical significance (p < 0.001) over static baselines, and Pearson correlation analysis ({rho} = 0.993 {+/-} 0.0002) confirms strict morphological fidelity. The full open-source codebase and benchmarks are publicly released. O_TBL View this table: org.highwire.dtl.DTLVardef@299dc5org.highwire.dtl.DTLVardef@3519a0org.highwire.dtl.DTLVardef@2586aborg.highwire.dtl.DTLVardef@1ac5610org.highwire.dtl.DTLVardef@1465c46_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 3:C_FLOATNO O_TABLECAPTIONQuantitative comparison on a common 60 s segment of Ninapro-like synthetic sEMG (single channel) with a 3 mV 50.3 Hz mains tone slightly drifted from the static notchs design centre at 50.0 Hz, stress-testing the adaptive corrector under a frequency mismatch. The Ninapro multi-channel aggregate (13.9 dB) reported in Section 3.4 uses mains exactly at 50 Hz (matched notch) and so achieves a higher {Delta} SNR. "MAC/sample" excludes the EMRMS square root and the pre-computed LMS sine/cosine. C_TABLECAPTION C_TBL

Abstract: Compartmental epidemic models conventionally treat the probability of moving between dis ease states as fixed over time, an assumption that sits uneasily with the reality of a pandemic in which lockdowns, mask mandates, vaccination roll-out, and the arrival of new variants continually reshape transmission. This paper develops a time-inhomogeneous Markov chain framework for the Susceptible-Exposed-Infectious-Removed (SEIR) process, in which each transition probability pab(t) is allowed to vary with calendar time while respecting the struc tural zeros implied by the SEIR compartmental flow. We derive the constrained maximum likelihood estimator of pab(t) under these structural constraints, establish its finite sample efficiency, asymptotic normality, and Wilson score confidence intervals, and construct a like lihood ratio test of the null hypothesis that a compartments exit probability is constant over time. We further propose a stochastic machine learning hybrid extension in which the raw, kernel smoothed transition probabilities are regressed on policy and mobility covariates using both a logistic generalized linear model and a random forest, allowing the framework to attribute time-inhomogeneity to observable interventions. The methodology is applied to a compiled daily, district level COVID-19 surveillance panel for Sierra Leone spanning March 2020 to December 2023 (16 districts, 1,401 days). The likelihood ratio test rejects time-homogeneity of the exposed to infectious transition in 15 of 16 districts and of the infectious-to-removed transition in 8 of 16 districts ( = 0.05), and the covariate augmented logistic model achieves an out of sample Brier score roughly 76 times smaller than a time homogeneous pooled baseline, with healthcare capacity and the time trend emerging as the most influential predictors in the random-forest component. These results provide statisti cal evidence that time-inhomogeneous, covariate informed Markov models offer a materially better description of district-level COVID-19 transmission in Sierra Leone than classical time homogeneous compartmental models, with implications for sub-national outbreak monitoring in resource constrained settings

Superspreading events (SSEs) produce extreme, rare bursts of disease transmission that standard compartment models, which assume population homogeneity, fail to capture. This inability to model heterogeneity in transmission rates can result in biased estimates of transmissivity. To address this limitation, we present a modular framework that treats SSEs as statistical outliers in case count time series and incorporates them into SIR-type models via pulse terms that transfer SSE cases directly from susceptible to infected compartments. This separation isolates anomalous SSE-driven transmission from background spread, which reduces bias when estimating mean transmission rates. We validate the approach on synthetic data generated by a stochastic model with embedded SSEs, demonstrating accurate recovery of the true non-SSE transmission parameter. We then apply the method to COVID-19 outbreaks in Hong Kong and the German district of Gutersloh, showing improved model fits and more robust estimates of background transmissivity both for a period with constant transmission and for a period with temporally structured NPI-driven heterogeneities. The framework's interchangeable outlier-detection, compartment, and SSE modules make it adaptable to diverse diseases and data contexts.

Approximate Bayesian computation sequential Monte Carlo (ABC-SMC) propagates its particles with a perturbation kernel, and with the standard Normal kernel it degrades sharply as the parameter dimension grows, a failure usually attributed to dimension itself. We show instead that it is governed by the quality of the summary statistics, with dimension entering only through a separate and milder mechanism, and that the two must act together for the Normal kernel to break. The first ingredient is covariance overinflation: the kernel covariance, estimated from the particle cloud, overshoots the true posterior covariance by a factor set by information loss in the summary statistics. We derive this overscaling factor in closed form for a Gaussian model with sufficient statistics and show that it stays modest at any dimension, shrinking toward its baseline value as the tolerance tightens; the extreme values seen in practice (of order 103) are a signature of insufficient summaries, not of dimension. The second ingredient is perturbation overconcentration: the normalised Normal step size concentrates around one as the dimension grows, so every proposal overshoots by the same factor. Either ingredient alone is harmless; only their combination breaks the Normal kernel. A Cauchy kernel (multivariate t with one degree of freedom) removes the concentration, keeping a positive acceptance rate under arbitrary overscaling at a bounded worst-case cost of 1.87x in expected squared jump distance. In a Metropolis-Hastings framework we derive closed-form acceptance rates for both kernels that illustrate the advantage of the Cauchy kernel in this limit. A series of full ABC-SMC computational experiments on five problems at d = 12, including a hierarchical gene-expression model, show the Cauchy reducing the sliced Wasserstein distance to the reference posterior by factors of up to 50 with the same simulation budget. Since the summary statistics are commonly insufficient for the models that require ABC, overinflation is structural and the Cauchy perturbation kernel is the right default for problems in higher dimensions.

Accurate estimation of the contrast transfer function (CTF) of tilt images is a critical first step in cryo electron tomography (cryoET), enabling reliable recovery of high-resolution structural information from thick, heterogeneous specimens. This challenge is especially acute in in situ cryoET, where macromolecules are imaged in their native cellular environment, often at high tilt and through substantial specimen thickness, with correspondingly low signal-to-noise ratios. Although CTF parameters can be later refined using reference-based approaches, accurate initial estimates are critical for downstream processing and the interpretability of tomographic reconstructions, yet they remain difficult to automate. Here, we present Quinoa, a software package designed to address these challenges. Quinoa first validates the tilt geometry and assesses data quality to generate robust initial estimates of defocus and phase shift. These estimates are then refined through optimization of a single global model, enabling precise fitting of the per-image defoci, tilt-dependent astigmatisms, time-dependent phase shifts, the specimen orientation (rotation, tilt and pitch) and the specimen thickness. Notably, and as a key distinguishing feature of this approach is that Quinoa fits equiphase-binned polar power spectra. This substantially reduces the computational cost of optimization without sacrificing accuracy, enabling more progressive and exhaustive refinement passes that further improve robustness. We validated Quinoa using both simulated and experimental data and benchmarked its performance against Warp, Ctfplotter, CTFMeasure, and AreTomo. Our results show that Quinoa is the most robust approach across all simulated cases, maintaining high accuracy even in the simultaneous presence of severe astigmatism, high specimen inclination and variable phase shift. Integrated recovery mechanisms further allow Quinoa to adapt automatically to a wide range of pixel sizes, defoci, astigmatisms and specimen thicknesses. Despite fitting a more complex and dynamic model, Quinoa remains extremely efficient due to extensive GPU acceleration, making it well suited for real-time monitoring during data collection as well as high-throughput offline batch processing. By improving automated CTF estimation in challenging tomographic data, Quinoa supports more accurate structural analysis of cells and tissues in situ.

Mechanistic ordinary differential equation models are widely used in systems biology to represent biochemical networks, population dynamics, cell-state transitions, and other biological processes; however, their predictive value depends critically on accurate parameter estimation from noisy and often sparse experimental data. In this tutorial, we present the Weak-form Estimation of Nonlinear Dynamics (WENDy) method as a forward-solver-free approach that reformulates parameter estimation as a covariance-corrected weak-form regression problem by integrating the model equations against compactly supported test functions. We present the background on the methodology through the lens of the familiar logistic equation, and we demonstrate applications of the method on real experimental data through two systems biology examples: a glycolytic oscillator with relatively dense time-course data and a sparse epithelial-mesenchymal cellstate transition model with multiple experimental replicates. Ultimately, using WENDy, we estimate interpretable biological parameters with uncertainty for systems with noisy and sometimes sparse available experimental data.

Significance: Continuous, noninvasive, and depth-resolved monitoring of blood-flow-related tissue dynamics remains an important unmet need. Speckle-contrast optical spectroscopy (SCOS), including interferometric implementations such as iSCOS, provides a scalable optical route to blood-flow sensing, but conventional continuous-wave approaches lack intrinsic depth selectivity. Time-of-flight (TOF) gating offers a way to separate superficial and deeper dynamic contributions in layered tissues, such as skin-muscle or scalp-cortex, by resolving photon path lengths. Aim: We introduce a swept-source, single-channel implementation of interferometric speckle-contrast optical spectroscopy (iSCOS) to obtain TOF-resolved temporal speckle contrast, {kappa}^2, from the measured field autocorrelation g_1, and evaluate its feasibility for depth-resolved blood-flow sensing. Approach: A swept-source iNIRS system operating at 780 nm acquired interferometric signals, which were Fourier-transformed along the optical-frequency axis to recover complex TOF-resolved speckle fields. Temporal speckle contrast was then estimated at each TOF gate indirectly from g_1 using the speckle-visibility relation. Diffusion-based numerical simulations were first used to compare the direct variance-based estimator and the indirect g_1-based estimator under varying reduced scattering coefficient, diffusion coefficient, additive noise level, and bi-layer geometry. Because the simulations showed that the g_1-derived {kappa}^2 estimator was substantially less sensitive to additive noise than the direct estimator, this estimator was used for the main phantom and in vivo analyses, while the direct estimator served as a simulation comparator. The g_1-derived estimator was then applied to liquid and bi-layer phantoms, followed by proof-of-concept in vivo measurements on the human forearm during cuff occlusion and on the forehead during a Sudoku task. Results: TOF-resolved kappa2 curves recovered with the g_1-derived estimator matched DWS theory across scattering coefficients, photon path lengths, and exposure times. The estimator preserved theoretical accuracy for additive noise amplitudes up to 50% of the field amplitude, whereas the direct variance estimator showed substantial noise-induced bias and required correction. Bi-layer simulations and phantom experiments reproduced the predicted direction and onset of TOF-dependent decorrelation-rate trends in layered media. In vivo, the recovered blood-flow index tracked the expected TOF-dependent cuff-occlusion and reactive-hyperemia response in the forearm. During the single-subject Sudoku task, the left-forehead recording showed a TOF-dependent relative blood-flow-index increase of +0.8 {+/-} 1.9% at TOF = 400 ps, +9.8 {+/-} 2.2% at TOF = 600 ps, and +15.2 {+/-} 5.6% at TOF = 800 ps. This pattern is consistent with increased sensitivity to deeper tissue at longer photon path lengths, but requires cohort-level validation before quantitative interpretation as cognitive activation. Conclusions: Coupling temporal speckle-contrast analysis with swept-source iNIRS yields a proof-of-concept, depth-resolved platform for blood-flow sensing. By estimating TOF-resolved speckle contrast through the g_1-derived {kappa}^2 route, TOF-iSCOS suppresses additive-noise bias while preserving sensitivity to deeper dynamic tissue layers. The present single-channel results bridge continuous-wave iSCOS, interferometric NIRS and time-domain diffuse correlation spectroscopy (TD-DCS), and motivate future multi-channel and cohort studies for scalable cortical hemodynamic monitoring.

Infectious disease transmission evolves as a dynamic process shaped by biological mechanisms, population behavior, and intervention policies, yet public health responses are often driven by lagging indicators. Accurate short- and long-term disease forecasting is essential for the timely deployment of intervention strategies, healthcare capacity planning, and uncertainty-aware, risk-informed decision-making. To address this challenge, three broad classes of forecasting models have traditionally been used: statistical, machine learning, and mechanistic approaches. However, each of these modeling paradigms faces fundamental limitations. In particular, traditional statistical models often lack the flexibility needed to capture complex disease dynamics, machine learning approaches require large, high-quality data streams, and mechanistic models are notoriously difficult to calibrate. To overcome these challenges, we propose a novel physics-informed machine learning (PIML) framework for forecasting infectious disease dynamics. Our approach simultaneously forecasts new case and hospitalization counts, along with other key epidemiological quantities such as the time-varying reproduction number. This is achieved through the design of a machine learning model and estimation strategy regularized by a system of differential equations that encode disease dynamics of the SIHR model, thereby bridging the gap between purely data-driven and mechanistic models. We demonstrate the proposed methodology through in-depth numerical studies and an application to COVID-19 data collected in the state of South Carolina.

Understanding how genomic surveillance capacity translates into population health outcomes is critical for designing effective pandemic response systems, yet the interaction between operational design and epidemiological dynamics remains insufficiently characterized. We develop an integrated analytical framework that links a whole-genome sequencing (WGS) - based surveillance network with a two - variant epidemiological transmission model to evaluate how surveillance operations influence variant detection, intervention timing, and health outcomes. The framework combines a modified susceptible - exposed - infectious - recovered - susceptible (SEIRS) model with a detailed operational representation of a centralized WGS surveillance network in India, incorporating sample collection, transport, batching, sequencing capacity, and reporting delays. We simulate 54 scenario combinations defined by three sequencing capacity levels, three sampling proportions, three variant emergence timings, and two variant profiles (high severity - high immune escape and low severity - low immune escape). Detection of a novel variant triggers a modeled intervention consisting of isolation of some diagnosed individuals, increased testing rates across disease states, and expanded access to hospitalization. Across simulations, the time from variant emergence to intervention implementation ranged from 73 to 351 days, depending on operational and epidemiological conditions. Increasing sampling proportion reduced detection time only when sequencing capacity was sufficient; under constrained capacity, higher sampling increased congestion and delayed detection. Expanding capacity from low to nominal levels substantially reduced turnaround times, with diminishing returns at higher capacity. Earlier detection consistently improved intervention effectiveness, with deaths averted ranging from 0.06% to 14.49% across scenarios. The cost per life - year saved ranged from INR 9,137 to INR 326,714 across all configurations, remaining below one to three times India ' s GDP per capita, consistent with established cost - effectiveness thresholds. These results demonstrate that the performance of genomic surveillance systems is jointly determined by operational and epidemiological dynamics. Effective surveillance design, therefore, requires coordinated optimization of sampling strategies and sequencing capacity to enable timely intervention and maximize population health benefits.

Under a range of behavioral and physiological conditions, spike times and local field potential (LFP) oscillations exhibit phase coupling within specific frequency bands. Classical measures such as spike-field coherence (SFC) and the phase-locking value (PLV) quantify this coupling but estimate the LFP spectrum independently of spike timing. We introduce Joint SSMT, a Bayesian state-space framework that jointly infers LFP spectrograms and spike-field coupling strength. The model treats narrowband LFP activity as a latent process evolving in continuous time, with spike trains linked to the complex spectral state through a Bernoulli-logistic model. In simulations, Joint SSMT accurately recovers coupling strength, denoises the spectrogram, and uses spike timing to resolve fine temporal structure in the LFP. Applied to propofol anesthesia data, the model identifies coupling at a specific slow-oscillation frequency where SFC and PLV report only broad low-frequency coupling. We extend Joint SSMT to trial-structured experiments and apply it to primate recordings during an associative learning task, revealing frequency-specific coupling in hippocampus and prefrontal cortex. We also derive closed-form expressions for SFC and PLV as functions of the generative model parameters. Across simulations and two primate datasets, Joint SSMT provides more frequency-specific coupling estimates with principled uncertainty quantification than classical PLV and SFC.

Respiratory diseases cause considerable morbidity in autumn and winter and are a priority in public health monitoring. In Germany, they are subject to a number of surveillance systems, including both pathogen-specific and syndromic indicators. In this paper we present a collaborative multi-target and multi-model real-time forecasting system rolled out during the 2024/25 season, and discuss differences to earlier efforts carried out during the COVID-19 pandemic. A total of nine models were run to generate forecasts of general practitioner consultations for acute respiratory infections (ARI), hospitalizations for severe acute respiratory infections (SARI) and confirmed cases of seasonal influenza and RSV. As all indicators were subject to retrospective revisions, forecasting models were combined with a nowcasting step. Whenever multiple models were available for the same indicator, we combined them into an ensemble. Nowcasts showed convincing performance, even though for some models Christmas break effects led to an upward bias in early January. Forecasts were overall well-calibrated and most models outperformed simple benchmark models. These improvements were generally more substantial for age-stratified than pooled targets, and concentrated at lead times of two to three weeks. Anticipating the peak timing and magnitude proved to be challenging, with many models predicting too flat curves with a too early turnaround (e.g. already in late January rather than mid-February for SARI). The combined ensemble forecast was among the best-performing approaches, but unlike in previous related projects did not consistently outperform individual models. We conclude by discussing learnings on the organization of collaborative forecasting projects in post-COVID-19 times and the potential of AI-supported modelling.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Microscopes are essential tools for discoveries on a scale invisible to the unaided human eye. The development of immuno-fluorescence followed by molecular biology techniques and fluorescent fusion proteins have revolutionised the use of optical microscopy in bioscience. The quality of the data produced is dependent upon the sample, its preparation and the instrument used. However, instruments can degrade over time without easily visible changes to the produced images and, in turn, negatively impacts results. By testing instruments and doing comparisons between results over time and between different instruments, problems can be highlighted and corrective action can be taken. Using small fluorescent beads the point spread function (PSF) of the microscope can be recorded and the image resolution measured. Beads were prepared in a concentration matched to the field of view size and dried onto coverslips and mounted on slides. The beads were then imaged as 3D Z-stacks of sufficient size to fully enclose the PSF of the system. This data was uploaded to OMERO and processed using OMERO-metrics, an OMERO plugin developed for this purpose. This paper summarizes the development of workflows and protocols to enable this process, presents the results obtained and demonstrates the detection of significant instrument issues.

The incidence rate ratio R is the standard measure for comparing event rates in clinical trials and epidemiology. In vaccine trials, the vaccine efficacy is VE = 1 - R. When events are rare, the two arm counts are Poisson. The estimator of R is heteroskedastic: its sampling variance changes with the data. So no fixed-width interval covers correctly everywhere. The usual log-Wald interval is undefined at zero events and covers poorly at small counts. Early vaccine and drug-safety readouts fall in exactly this regime. We show that a single reparameterization collapses this bivariate problem to an effective one-parameter family with a quadratic variance function, whose variance-stabilizing transformation is 2 arcsinh(sqrt(R)). The reduction yields a closed-form confidence interval for R. Its two leading errors, a curvature bias and the variability of the estimated scale, each admit a closed-form correction with no tuning constants. In a Monte Carlo study of our seven arcsinh variants and five competitors, the +Curve+Stu variant covers within 0.002 of the nominal 0.95 for about 50 control and 5 treatment events. Its width is on par with the best competitor. It avoids the conservatism and zero-count breakdown of log-Wald and MOVER. For moderate counts, we recommend this interval; for sparser data, our Bar-Lev and Enis count-shift variant is more robust. The result is a ready-to-use, closed-form interval for the low-count regime. We illustrate it on early Covid-19 vaccine-efficacy readouts and provide reference implementations in R and Python.

The measurement of light received at the cornea of the eye is a paramount consideration for the understanding of the relation between environmental illumination and the non-image-forming effects of light. The field of view (FOV) at the cornea is less than a full hemisphere, because it is partially occluded by human facial morphology. The International Commission on Illumination (CIE) has defined a standard model of human FOV. A suitably designed physical occluder attached to the sensor (of a light meter) has been proposed as a means of incorporating the effect of human FOV when taking measurements. Similarly, when using simulation to predict light received at the cornea, a geometrical description of the occluder at the eye point(s) can be added to the 3D model of the scene. The first occluder model proposed to represent CIE human FOV was enumerated in terms of: the CIE definition; the radius of the occluder; and, the radius of the light sensor disc. We present a simpler model based only on the CIE definition and the occluder radius. Both models were tested using a virtual goniophotometer. Various sensor response functions describing the spatial sensitivity across the sensor disc, including several we characterized through laboratory measurements, were included in the test. For all functions considered, the performance of the simpler occluder model was equivalent to or better than the model first proposed.

AO_SCPLOWBSTRACTC_SCPLOWEarly tumor detection is critical for improving patient survival and recovery. Clinically, tissue palpation is routinely used to identify regions of abnormal stiffness, a hallmark of many pathological conditions. However, palpation is restricted to anatomically accessible sites and remains highly operator dependent. Here, we introduce a method for real-time quantitative stiffness mapping using an unmodified commercial endoscope, with the goal of enhancing diagnostic capabilities and restoring mechanical feedback during endoscopic procedures. Our approach combines shear wave elastography with speckle imaging and an innovative synchronization strategy that enables the measurement of shear wave propagation using an unmodified commercial endoscope. The resulting wave fields are analyzed with the noise-correlation-inspired (NCI) method[1], providing pixel-wise estimates of shear wave velocity and, consequently, quantitative maps of local tissue stiffness. The method demonstrated robust performance in both benchtop and endoscopic configurations. Validation was achieved on polymer phantoms as well as on ex vivo and in vivo biological tissues, highlighting its potential for minimally invasive biomechanical imaging and real-time tissue characterization.

Flow imaging microscopy (FIM) is an important technology for high-throughput characterization of microscopic particles and microorganisms. However, conventional FIM relies on single-plane imaging (SPI), resulting in out-of-focus particles, reduced measurement precision, and incomplete characterization of irregularly shaped objects extending along the z-axis. To address these limitations, a volumetric flow imaging (VFI) framework was developed and implemented on the portable ARTiMiS platform. This approach captures multiple frames along the z-axis and extracts the highest fidelity image for each particle, which can also be used for single image generation with all particles in focus (i.e., all in focus image) and for three-dimensional reconstruction of irregularly shaped objects. Benchmarking VFI with microspheres, live cells (Chlorella vulgaris), and filamentous cyanobacteria demonstrated increased fraction of particles in focus, reduced variability in particle size measurement, and increased resolvability of elongated particles in comparison to conventional SPI on commercially available FIM technologies. For C. vulgaris, VFI-derived size distributions closely matched curated FlowCam measurements without requiring post-processing to exclude out-of-focus particles. All-in-focus image reconstruction enabled simultaneous visualization of particles distributed across multiple depths and consistently resolved a greater proportion of filamentous structures as compared to SPI. For Aphanizomenon sp., Dolichospermum sp., and Planktothrix agardhii, the SPI approach captured only 84%, 61%, and 58%, respectively, of the total filament length resolved by AIF reconstruction. Beyond image-based characterization, VFI enabled estimation of dynamic particle properties such as sinking velocity and mass density. Application of this framework to C. vulgaris cultures revealed distinct mass-density trajectories under nitrogen-replete and nitrogen-deplete conditions, with cell mass density increasing over time under nitrogen-replete conditions and decreasing under nitrogen deprivation. Collectively, these results establish VFI as a next-generation framework for FIM that expands its analytical capabilities beyond conventional morphometric characterization and provides new opportunities for single-cell-enabled environmental monitoring and biomanufacturing.

Optogenetics enables experiments in out-of-equilibrium conditions to clarify biological mechanisms and quantify biophysical parameters. However, modelling and control techniques to study mammalian cell biology under optogenetic perturbation remain underutilised. Here, we benchmark these methods within mammalian cells by steering nucleocytoplasmic transport via the optogenetic LEXY protein in outcome-driven microscopy. First, we employ system identification to obtain models that predict transport dynamics by minimising the prediction error. We quantify this prediction accuracy for one biophysical model and two black-box models. Second, we evaluate closed-loop control efficacy by steering transport along a predefined trajectory using model-free Proportional Integral (PI) control, model-based Linear Quadratic Regulation (LQR) and Model Predictive Control (MPC). Both the predictive models and the applied control techniques demonstrate robust performance against cell-to-cell variation. This biological variation is quantified by the parameter distributions obtained from model identification with single-cell trajectories. While we show that model-free techniques such as PI and gain-scheduled PI achieve steering without explict model knowledge, predictive architectures offer better performance under this cell-to-cell variation and time-varying setpoints. Moreover, black-box predictive accuracy suggests that this model-based control is possible, even when explicit mechanistic understanding is missing. Ultimately, we demonstrate that predictive modelling and optogenetics enable quantitative characterisation and precise manipulation of mammalian cells, while offering practical guidelines for the implementation of these techniques.

Cardiac blood flow is a regulator of several important developmental and remodelling processes in the heart, including through fluid shear forces sensed by the endothelial cells lining the heart. However, optically mapping these flow fields in the complex 3D geometry of the heart is challenging even in transparent animal models such as the zebrafish. One of the main challenges is the difficulty in measuring the out-of-plane (axial) velocity component, preventing accurate mapping of the complete 3-component-3-dimension (3C-3D) blood flow velocity field; image-based techniques such as microscopic particle image velocimetry ({micro}PIV) traditionally only provide the in-plane flow components. Here we present a computational approach to achieve full time-varying 3C-3D blood flow vector mapping using a standard selective plane illumination microscope (SPIM), based on robust cardiac phase assignment, precise measurement-driven registration of sequentially acquired z-stacks, and PIV data fusion from multiple sample orientations. Our approach holds the key to understanding the complex dynamic flow fields within the developing heart, and their role in shaping cardiac development.